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I. INTRODUCTION 

Over the last few years, many non-trivial quantum phenomena have been discovered in the low- 
dimensional magnetic materials. One of the most actively studied problems in the condensed matter 
physics at present is a quantum tunneling of magnetization, i. e. the effect, which has been recently 
developed in a large number of single- molecule magnets (see Ref. [1] and references therein). By the term 
single-molecule magnets, one denotes the small clusters of magnetic metal ions that usually possess an 
extraordinary strong magnetic anisotropy. Hence, the single-molecule magnets often provide very good 
examples of so-called Ising-like spin systems with a strong uniaxial magnetic anisotropy. Of course, the 
Ising anisotropy by itself cannot be a source of the quantum spin tunneling experimentally observed in 
these systems. It turns out, however, that this quantum phenomenon arises in the most cases due to 
the higher-order crystal-field terms. According to a number of experimental and theoretical studies it 
is now quite well established, that the observed spin tunneling originates to a major extent from the 
second-order biaxial crystal-field potential, at least in Fe4^, Feg^, Fcig^, or Mn4^ compounds. 

The immense interest in the magnetic properties of small magnetic clusters shed light on the effect of 
single-ion anisotropy terms D (uniaxial anisotropy) and E (biaxial, also called rhombic anisotropy). In 
contrast to the quite well understood role of the both single-ion anisotropics D and E in the small mag- 
netic clusters (zero-dimensional systems) , the situation is much more complicated and also obscure in one- 
and two-dimensional spin systems. In fact, the ground-state properties of a spin-S* Ising model with the 
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rhombic crystal-field potential E, have been only recently examined by Oitmaa and von Brasch within an 
effective mapping to the transverse Ising model^' . On the basis of this effective mapping, the T = quan- 
tum critical point can be exactly located for the one-dimensional model, while for the two-dimensional 
models they can be obtained with a high numerical accuracy using the linked-cluster expansion method^'''. 
Nevertheless, the finite temperature behaviour of these models has not been investigated in detail beyond 
the standard mean- field and effective- field theories^, random phase approximation^, or linked cluster 
expansion^°. It should be stressed that the biaxial anisotropy essentially influences the magnetic prop- 
erties of a large number of polymeric molecular-based magnetic materials, too. Prom the most obvious 
examples one could mention: NiFa", NiNOs-eHaO^^, Ni(CH3COO)2.4H20i3, Mn(CH3COO)2.3H20", 
CoF2^^, CoCl2.6H20^^ and a series of compounds Fe(dc)2X^^, where X stands for halids and dc for the 
dithiocarbamate or diselenocarbamate groups, respectively. 

Owing to this fact, in this article we will focus on the uniaxial and biaxial crystal- field anisotropics 
affecting the magnetic behaviour of the mixed spin-1/2 and spin-1 honeycomb lattice. By assuming an 
Ising-type exchange interaction between the nearest-neighbouring spins, the model under investigation can 
be exactly treated through the generalized star-triangle mapping transformation. The considered model 
thus provides a noble example of the statistical system, which enables to study an interplay between 
quantum effects and temperature in a spontaneously ordered magnetic system. Moreover, the magnetic 
structure of a mixed-spin honeycomb lattice occurs rather frequently also in the molecTilar magnetism, 
what clearly demonstrates a large family of polymeric two-dimensional compounds of chemical formula: 
A^M"m"^(C204)3^^, where stands for a non-magnetic univalent cation N(C„H2„+i)4 or P(C„H2„+i)4 
(n = 3 - 5), M" and M"^ denote two- and threc-valent metal atoms Cu"(5' = 1/2), Ni"(S' = 1), 
Co"(5' = 3/2), Fe"(S' = 2) or Mn"(5' = 5/2) and Cr"i(S' = 3/2) or Fc^^^S = 5/2), respectively. Indeed, 
the crystal structure of these polymeric molecular-based magnetic materials consists of the well-separated 
two-dimensional layers in which regularly alternating M^^ and M^^^ magnetic metal atoms constitute more 
or less regular honeycomb lattice (Fig. 1). As a consequence of the anisotropic crystalline structure of 
these materials, one should also expect a relatively strong uniaxial (Ising-like) anisotropy, as it has 
already been suggested in the theoretical studies based on the effective-field theory and Monte-Carlo 
simulations^^. Hence, the magnetic compounds from the family of oxalates A^M"M^^^(C204)3 represent 
good candidates to be described by the proposed model. 

The outline of this paper is as follows. In the next section the detailed description of the model system 
will be presented and then, some basic aspects of the transformation method will be shown. Section 3 
deals with the physical interpretation of the most interesting results and finally, some concluding remarks 
are drawn in Section 4. 

II. MODEL AND METHOD 

Let us consider the magnetic structure of a mixed-spin honeycomb lattice schematically depicted in 
Fig. 1. To ensure exact solvability of the model under investigation, we will further suppose that the sites 
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of sublatticc A arc occupied by the spin-1/2 atoms (depicted as full circles), in contrast to the sites of 
sublattice B that are occupied by the spin-1 atoms (open circles). By assuming the Ising-type exchange 
interaction J between nearest-neighbouring spins, the total Hamiltonian of the system takes the following 
form: 

3JV JV JV 

n = jj2 sr^'j + DY^^sif + EY,[{stf - [sif], (1) 

(fe,i) fees fees 

where iV is a total number of sites at each sublattice, /t| and S'^{a = x,y,z) denote the standard 
spatial components of the spin-1/2 and spin-1 operators, respectively. The first summation in Eq. (1) is 
carried out over the nearest-neighbouring spin pairs, while the other two summations run over the sites 
of sublattice B. Apparently, the last two terms D and E are the crystal-field potentials that measure a 
strength of the uniaxial and biaxial anisotropy acting on the spin-1 atoms. It is also worth noticing that 
there is one-to-one correspondence between the Hamiltonian (1) and the effective spin Hamiltonian with 
three different single-ion anisotropies D^, and D^: 

3N N N N 

H = JY1 -^IAI + ^(5|)2 + 5^(4^)2 + Dy J2iSir- (2) 

(fe,j) fees fees fees 

In fact, one can easily prove the equivalence between the two effective spin Hamiltonians (up to the 
unimportant additive constant, for a comparison see^*^), which can be achieved using this simple mapping 
between the relevant parameters included in the Hamiltonians (1) and (2), respectively: 

D = D' -^{D'' +Dy), and E=^{D=' - D^). (3) 

It should be also mentioned here that by neglecting the biaxial anisotropy, i. e. setting £" = in Eq. 
(1) or equivalently = in Eq. (2), our model reduces to the exactly soluble model of Goncalves^"'^. 
Accordingly, in this work we will in particular examine the effect of biaxial anisotropy on the thermody- 
namical and dynamical properties of the model imder consideration. Nevertheless, the E term emerging 
in the Hamiltonian (1) should cause non-trivial quantum effects, since it introduces the x and y compo- 
nents of spin operators into the Hamiltonian and thus, it is responsible for the onset of local quantum 
fluctuations that are obviously missing in the Ising model with the uniaxial crystal-field potential D only. 

It is therefore of interest to discuss the origin of biaxial anisotropy. The origin of this anisotropy term 
consists in the low-symmetry crystal field of ligands from the local neighbourhood of spin-1 atoms. A 
threefold symmetry axis oriented perpendicular to the honeycomb layer, however, prevents the appear- 
ance of biaxial crystal-field potential in a regular honeycomb lattice with a perfect arrangement of the 
oxalato groups, as well as magnetic metal atoms. On the other hand, the small lattice distortion, which 
occurs rather frequently in the low- dimensional polymeric compounds due to the Jahn- Teller effect, can 
potentially lower the local symmetry. In consequence of that, the distortion of lattice parameters can 
be regarded as a possible source of the biaxial anisotropy. The most obvious example, where the lattice 
distortion removes the threefold symmetry axis represents the single- molecule magnet Fe4, in which three 
outer Fe atoms occupy two non-equivalent positions around one central Fe atom^. 
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Let us turn our attention to the main points of the transformation method, which enables an exact 
treatment of the model under investigation. Firstly, it is very convenient to write the total Hamiltonian 
(1) as a sum of the site Hamiltonians Tik'- 

N 

W=E^'^' (4) 

where each site Hamiltonian Hk involves all interaction terms associated with the appropriate spin-1 
atom residing on the fcth site of sublattice B: 



Hk = SIE, + [SlfD + [{Slf - {Slf]E, 



(5) 



with E}~ = J (fill + Afe2 i^ks)- While the Hamiltonians (5) at different sites commute with each other 
{[HijHj] — 0, for each i ^ j), the partition function of the system can be partially factorized and 
consequently, rewritten in the form: 

JV 

Z = Tr^^y II TVS, exp(-/JWfc). (6) 

k=l 

In above, (3 = l/{kBT), ks is Boltzmann's constant, T the absolute temperature, Tr|^} means a trace 
over the spin degrees of freedom of sublattice A and Trg^ stands for a trace over the spin states of fcth 
spin from sublattice B. So, a crucial step in our procedure represents the calculation of the expres- 
sion Trg^ exp(— With regard to this, let us write the site Hamiltonian (5) in an usual matrix 
representation: 

f D + Ek E \ 

(7) 



V 




E 




D-Ek J 

in a standard basis of functions | ± 1), |0) corresponding, respectively, to the three possible spin states 
5^ = ±1, of fcth atom from sublattice B. Obviously, it is easy to find eigenvalues of the site Hamiltonian 
(7), however, with respect to further calculation, it is more favourable to obtain directly the matrix 
elements of the expression exp(— /JW/j). Using the well-known Cauchy integral formula, one readily 
obtains the matrix elements for an arbitrary exponential function of the site Hamiltonian (7): 

/ cosh(ae) + %- sinh(ae) | sinh(ae) \ 



exp{a'Hk) = exp(a£>) 



V 





f sinh(ae) 



ex.p{—aD) 
cosh(ae) ^ 



(8) 



Q- sinh(ae) / 



where O = \JE^-\- E"^ and a marks any multiplicative function. After substituting a — —j3 into the 
Eq. (8), the calculation of the relevant trace Tr^^ exp(— /JTife) can be accomplished, moreover, its explicit 
form immediately implies a possibility of performing a standard star-triangle mapping transformation: 



TVS, exp(-/3Hfc) = 1 + 2 exp(-/3Z)) cosh(/3y'j2(/,|, + + M^s)' + 



= A exp 



/?^(MfclMfe2 + Mfe2Mfe3 + Mfe3Mfel) 



(9) 
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which replaces the partition function of a star, i. c. the four-spin cluster consisting of one central spin-1 
atom and its three nearest-neighbouring spin-1/2 atoms, by the partition hmction of a triangle, i. e. 
the three-spin cluster comprising of three spin-1/2 atoms in the corners of equilateral triangle (see Fig. 
1). The physical meaning of the mapping (9) is to remove all interaction parameters associated with 
the central spin-1 atom and to replace them by an effective interaction R between the outer spin-1/2 
atoms. It is noteworthy, that the both mapping parameters A and R axe " self-consistently" given by the 
transformation equation (9), which must be valid for any combination of spin states of three spin-1/2 
atoms. In consequence of that one obtains: 

A=(^^,<^i)"\ /3i?=ln(|i), (10) 

where we have introduced the functions $i and $2 to write the transformation parameters (10) in more 
abbreviated and elegant form: 



$1 = 1 + 2 exp(-/3£)) cosh (^(3^ {3J/2y + E^) , 



$2 = 1 + 2 exp(-/3D) cosh (3 J ( J/2)2 + E^]. (11) 



When the mapping (9) is performed at each site of the sublatticc B, the original mixed-spin honeycomb 
lattice is mapped onto the spin-1/2 triangular lattice with the effective interaction R given by the "self- 
consistency" condition (lO)-(ll). As a matter of fact, the substitution of the mapping transformation (9) 
into the partition function (6) establishes the relationship: 

Z{(3,J,D,E) = A^Zt{(3,R), (12) 

between the partition function Z of the mixed-spin honeycomb lattice and the partition function Z^ of the 
corresponding spin-1/2 triangular lattice. Above equation constitutes the basic result of our calculation, 
since it enables relatively simple derivation of all required quantities such as magnetization, quadrupolar 
moment, correlation function, internal energy, specific heat, etc. Moreover, by combining (12) with (9) 
one easily proves the validity of following exact spin identities: 

(/i(Mf,M|,-,MD) = (/i(/«i>M|,-,Mfe))t, 

lf(Qx qV qz ,,z ,,z ,z \\ / Trg^,/2(6'g,6'fc,S'fc,/ifcl,/Xfc2,Mfc3)6Xp(-/3'Hfc) \ , , 

\ Trs^ exp(-/3Wfe) / 

where (...) represents the standard canonical average over the ensemble defined by the Hamiltonian (1) 
and {...)t canonical average performed on the spin-1/2 Ising triangular lattice with the effective exchange 
interaction R (lO)-(ll). Furthermore, /i is an arbitrary function of the spin variables belonging to 

the sublattice A, while /2 denotes an arbitrary function depending on the fcth spin from sublattice 
B and its three nearest-neighbours from sublattice A. Applying the first of spin identities (13), one 
straightforwardly attains the following results: 

mA = {fiki) = {f^ki)t=mt, (14) 
CA = (AfeiAfe2) = (AfciAfc2)t = ct, (15) 

tA = (Afel Afe2Afe3) = (AfclAfc2Afc3)t = tt, (16) 
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while the second of spin identities (13) enables a derivation of quantities depending on the spin variable 
from sublattice B, as well: 

mB = {Sl) = -3mA{K, + K2)/2 - 2tA{K, - 3K2), (17) 

9l = ((^fe)') = {K5 + 3i^6)/4 + Sca{K5 - Ke), (18) 

ql = {{Slf) = {Kr + 3Ks)/4 + 3cA{Kr - Ks), (19) 

Q% = {{SD') = {K3 + 3Ki)/A + 3cA{K3-Ki). (20) 

In above, ttia {m,B) labels the single-site magnetization at sublattice A (B), q'^{a = x,y,z) are different 
spatial components of quadrupolar moment and finally, ca and Ia static pair and triplet correlation 
functions between the relevant spins of sublattice A, respectively. Obviously, an exact solution for the both 
sublattice magnetization and quadrupolar moment require the knowledge of the single-site magnetization 
TOj, nearest-neighbour pair correlation hmction Cj and triplet correlation hmction tt on the corresponding 
spin-1/2 triangular lattice unambiguously given by (lO)-(ll). Fortunately, the exact solution for these 
quantities on spin-1/2 triangular lattice are known long time ago, hence, one can utilize the final results 
from references'^. Finally, the coefScients emerging in the previous set of Eqs. (17)-(20) are listed below: 

K,=F,{3J/2), K2 = Fi{J/2), i^a = i^2(3J/2), K^ = F2{J/2), 

K5=F3{3J/2,-E), Ke = Fs{J/2,-E), Kr = F3{3J/2,E), Ks = Fs{J/2, E), (21) 
where we have defined the functions Fi{x), F2{x) and F3{x,y) as follows: 



^ , , X 2sinh(/3Va;2 -h£;2) 

Fi{x)- 



F2{x)-- 



' ^x^ + E'^ exp(/3£)) + 2cosh(/3\/a;' + i;^) ' 

2cosh(/3Va;^ + j;^) 
' exp(/3D) 2 cosh(/3yP+l^' 



exp(/3D) + cosh(/3v/a?2~fy2) y sinh(/3^a;2 -f 

r'j,\x,y) = , — I , , 

exp(/3i:') + 2 cosh(/3 V^a;2 + j/2) + yi cxp{(3D) + 2 cosh(/3A/a;2 + y^) 

At the end of this section, we will also provide an exact result for one dynamical quantity - time- 
dependent autocorrelation function. It should be noted here that the exactly soluble models offer only 
seldom the possibility to investigate their spin dynamics. On the other hand, the dynamical quantities 
such as autocorrelation and correlation functions are important also from the experimental point of view, 
because their magnitude directly determines the scattering cross section measured in the inelastic neutron 
scattering experiments^'', or the spin-lattice relaxation rate provided by the nuclear magnetic resonance 
(NMR) techniques'^. In this work, an exact treatment for the time-dependent autocorrelation function 
will be elaborated. As a starting point for the calculation of the autocorrelation function C^^^g{t) can 
for convenience serve the second of exact spin identities (13): 

2\ Trs, exp(-/3H,) /' 
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where the symmetrized form in the definition of C^^j^ is used to construct a hermitian operator, 
S^{t) — cxp( cxp(— ) represents the Heisenberg picture for the time-dependent operator 
S^{t), h stands for Planck's constant and i = \/"-T. Next, the matrix representation of exp(± ) 
can be readily obtained by putting a = ±^ into Eq. (8). Then, after straightforward but a little bit 
tedious calculation, one arrives to the final result for the dynamical autocorrelation function: 



+^H4 - ^^*) oWTe^ - 



III. RESULTS AND DISCUSSION 



Before proceeding to the discussion of the most interesting results, it is noteworthy, that the results 
derived in the previous section are rather general, i. e. they are valid for the ferromagnetic (J < 0), 
as well as ferrimagnetic (J > 0) version of the model under consideration. In what follows, we will 
restrict ourselves to the analysis of the ferrimagnetic model only, since the polymeric compounds from the 
family of oxalates^^ fall mostly into the class of ferrimagnets. Nevertheless, it appears worthwhile to say 
that magnetic behaviour of the ferrimagnetic system completely resembles that one of the ferromagnetic 
system. Finally, it should be emphasized that the mapping (9) remains invariant under the transformation 
E ^ —E. As a result, one may consider without loss of generality the parameter > and consequently, 
X-, y- and 2;-axis then represent the hard-, medium- and easy-axis for a given system. 



A. Ground-state properties 



At first, we will take a closer look at the ground-state behaviour. Taking into account the zero- 
temperature limit T — » 0+, one finds following condition for a first-order phase transition line separating 
the magnetically ordered phase (OP) from the disordered phase (DP): 



From Eqs. (16)-(22), moreover, one easily attains analytical results for the single-site sublattice magne- 
tization {tha, tub), total single-site magnetization m = (m^ +mB)/'i and different spatial components 
of the quadrupolar moment q%{a = x, y, z) in the both phases, as well: 

1 ^ 1 2 

OP: mA = — -, ruB = — , ^ , m = — ^ 



2 /7\2 / X 2 ' 4 



= \il- , ' ), Ql = \i^+ I ' ). al = 1.0; (26) 
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DP: mA = 0.0, ms = 0.0, to = 0.0, 

g| = 1.0, ql = 1.0, g^ = 0.0. (27) 

For better illustration, Fig. 2 depicts the ground-state phase diagram in the E/J-D/J plane (Fig. 2a) 
together with the zero-temperature variations of the magnetization and quadrupolar moment in the OP 
(Fig. 2b, the value of uniaxial anisotropy D/J = 0.0 has been chosen not to pass through the phase 
boundary). It is worthy to mention that by neglecting the biaxial anisotropy, i. e. setting E/J = 0.0, 
one recovers from the phase boundary condition (25) a boundary uniaxial anisotropy D/J = 1.5, which 
has been already reported by Gongalves several years ago^^. In this limit, the OP corresponds to the 
simple ferrimagnetic phase in which both sublattice magnetization are fully saturated and also antiparallel 
oriented with respect to each other (in fact, niA = —0.5 and thb = 1-0). 

The situation becomes much more complicated by turning on the biaxial anisotropy E. Even though the 
sublattice magnetization niA remains at its saturation value in the whole OP, the sublattice magnetization 
niB is gradually suppressed by increasing the biaxial anisotropy strength. In contrast, neither sublattice 
magnetization, nor the quadrupolar moment do not depend within either ground state phase on the 
uniaxial crystal-field potential D. Of course, the relevant change of sublattice magnetization uib must 
reflect a violation of a perfect ferrimagnetic spin arrangement in the OP. To achieve the non-saturated m b 
at T = 0, some spins of sublattice B must flip from the | -f- 1) to | — 1) and/or |0) state(s). It is therefore 
of great importance to identify the magnitude of the quadrupolar moment qg. Since the quadrupolar 
moment approaches in the OP its saturation value = 1.0 independently oi E/J, a presence of the |0) 
states can be thus clearly excluded. These observations would suggest, that the biaxial anisotropy causes 
in the OP a spin tunneling between the | -|- 1) and | — 1) states, whereas the stronger the ratio E/J, the 
greater the population of the | — 1) state. Anyway, the probabilities to find the spin-1 atom in the | ± 1) 
state are given by these simple expressions: p(| ± 1)) = . Altogether, the spin configuration referring 

to the OP at T = can be characterized as follows: all spin- 1/2 atoms are wholly ordered in their spin 
down positions (m^ = —0.5), while the spin-1 atoms occupy with the probability p(|±l)) either the | + 1), 
or I — 1) state. It should be also pointed out, that the condition q^ > q'^ is always satisfied when E > 0. 
This inequality between the spatial components of quadrupolar momentum provides a confirmation, that 
X- and 2/-axis represent under the assumption E > the hard- and medium-axis in the OP. 

At last, let us consider the spin ordering within the DP. Interestingly, the DP remains unaltered 
no matter whether the biaxial anisotropy is zero, or not. Indeed, all spin-1 atoms occupy in the DP 
exclusively the |0) state, because of ttib = qs = O-O- Contrary to this, the components of quadrupolar 
moment perpendicular to the z-ssds acquire in the DP their maximum value q^ = q^B ^ ^-^^ These results 
can be thought as an independent check for the scenario that accompanies the phase transition from the 
OP to DP: all spin-1 atoms indeed tending to align into the x — y plane. Accordingly, the magnetic 
order is completely destroyed, in fact, the vanishing magnetization niA implies a state of complete spin 
randomization at sublattice A. Therefore, the DP does not exhibit any long-range magnetic order even 
at T = 0. 
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Now, another interesting question arises, namely, whether the spin-1 atoms can fluctuate in the OP 
between their allowable \ ± 1) states. In order to obtain a reliable answer to this question, the time- 
dependent autocorrelation function (24) will be analysed. In the zero-temperature limit, the dynamical 
autocorrelation function C^^^q gains after straightforward calculation: 

(|)^ + (f) ^cos(^y^( |)^ + (f)^) 

(i)^ + (f) 

which in turn proves that C^^j^ is periodic in time with the angular frequency a;„ = x^/ (f)^ (t)^ 
the recurrence time r = — ^====. According to Eq. (28), the dynamical autocorrelation function 
does not depend in the ground state on the uniaxial anisotropy D. Owing to this fact, we will further 
neglect this anisotropy parameter and set D/J = 0.0. For illustrative purposes, the time variation of 
the autocorrelation function C^^j^ are displayed in Fig. 3 for several values of the biaxial anisotropy 
E/J = 0.1, 0.5, 1.0 and 2.0. It appears worthwhile to make a few remarks on foregoing results. Since 
the autocorrelation function varies in time, it clearly demonstrates the zero-temperature spin dynamics 
between the allowable \ ± 1) states. From the analytical solution (28) as well as depicted behaviour one 
can moreover deduce a physical interpretation of the spin dynamics, namely, the spin system necessarily 
recovers after the recurrence time t always its initial state, whereas the stronger the ratio E/J, the shorter 
the recurrence time r. In addition, the increasing strength of the biaxial anisotropy enhances also the 
time- variation of C^^j^ (i. e. the amplitude of oscillation). This result is taken to mean, that increasing 
biaxial anisotropy enlarges also a number of the spin-1 atoms, which tunnel during the recurrence time 
between the | ± 1) states. Since the equilibrium magnetization docs not varies in time, a number of atoms 
that tunnel from | + 1) to | — 1) state, must be definitely the same as a number of atoms that tunnel 
from the | — 1) to | -|- 1) state. These findings have an obvious relevance to the understanding of the 
zero-temperature spin dynamics, because they enable its explanation from the microscopic viewpoint. 



B. Finite-temperature behaviour 



In this part, we would like to make some comments on the finite- temperature behaviour of the system 
under investigation. Let us begin by considering the effect of uniaxial and biaxial anisotropies on the 
critical behaviour. For this purpose, two typical finite-temperature phase diagrams are illustrated in 
Fig. 4a and 4b. In both figures, the OP can be located below the phase boundaries depicted as solid 
lines, while above the relevant phase boundaries the usual paramagnetic phase becomes stable. A closer 
mathematical analysis reveals, that the temperature-driven phase transition between these two phases 
is of second-order and belongs to the standard Ising universality class. More specifically. Fig. 4a shows 
the critical temperature as a function of the uniaxial anisotropy D/J for several values of the biaxial 
anisotropy E/J. The dependence critical temperature versus uniaxial anisotropy is quite obvious, when 
increasing D/J, the critical temperature tends monotonically to zero as many as the boundary value (25) 
is achieved. While the anisotropy term D forces the spins to lie within x — y plane when D > 0, the E 
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term tries to align them into y — z plane. Accordingly, the increasing strength of the biaxial anisotropy 
supports the magnetic long-range order related to the OP when D / J > 1.5 and hence, it survives mitil 
stronger anisotropics D/ J. As far as the region Z? < is concerned, the biaxial anisotropy substantially 
lowers the critical temperature of the OP. Apparently, this behaviour arises as a consequence of the 
fact, that the E term simplifies the transition between the | ± 1) states due to the non-zero quantum 
fluctuations. Thus, one can conclude that the quantum fluctuations macroscopically manifest themselves 
in the reduction of the critical temperature for the easy-axis uniaxial anisotropy (i. e. for D <0, where 
the model Hamiltonian (1) works extremely well). 

To illustrate the influence of the biaxial anisotropy on the critical behaviour, the critical temperature 
versus biaxial anisotropy dependence is shown in Fig. 4b for several values of the uniaxial anisotropy. 
As one would expect, the critical temperature gradually decreases with increasing the biaxial anisotropy 
strength for any D <0. In agreement with the aforementioned arguments, the appropriate depression of 
the critical temperature can be again attributed to the quantum fluctuations, which become the stronger, 
the greater the ratio E/J. Apart from this rather trivial finding, one also observes here the interesting 
dependences with the non-monotonical behaviour of the critical temperature. Namely, for D/J » 0.0 
the critical temperature firstly increases and only then gradually decreases with the biaxial anisotropy 
strength (see for instance the curve for D/J = 1.3). To explain such a behaviour, it should be realized 
that the spin-1 atoms are preferably thermally excited to the |0) state when D > 0, what means, that they 
arc preferably excited to the x — y plane. Since the biaxial anisotropy tries to align them into the y — z 
plane, it favors the long-range order along 2:-axis in that it prefers the spin tunneling between the | ± 1) 
states before the population of the |0) one. The most interesting result to emerge here is that there is a 
strong evidence, that aforementioned argument explains an existence of the OP even under assimiption 
of extraordinary strong anisotropics D/J > 1.5. In fact, the magnetic long-range order related to the OP 
occurs under this condition for the strong enough biaxial anisotropics only. Surprisingly, the magnetic 
long-range order results in such a peculiar case from the quantum fluctuations (spin tunneling) caused 
by the biaxial anisotropy. 

Now, let us provide an independent check of the critical behaviour by studying the thermal dependences 
of magnetization. The single-site magnetization against the temperature are plotted in Fig. 5 for the 
uniaxial anisotropy D/J = —2.0 and several values of the biaxial anisotropy E/J. Fig. 5a shows a typical 
situation observed by turning on the biaxial anisotropy E/J: the greater this anisotropy parameter, the 
stronger the reduction of sublattice magnetization ms due to the | ± 1) spin tunneling. As it is apparent 
from this figure, the total magnetization exhibits in general the standard Q-type dependences. The most 
striking thermal variations of the total magnetization can be evidently found for the biaxial anisotropies 
close to the value E^/J = \/27/2, at which niA fully compensates ms in the ground state (Fig. 5b, see 
also Fig. 2b). It turns out, however, that all marvellous thermal dependences of the total magnetization 
stemming from the identical origin - the magnetization of sublattice A is thermally more easily disturbed 
than the magnetization of sublattice B. As a result, the P-type dependences of total magnetization occur 
for E < E^, when the prevailing magnetization ms exhibits smaller thermal variation than the lower 
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magnetization rriA (sec cases E/J = 2.58 and 2.59 in Fig. 5b). Based on our earlier remark concerning 
the groimd-state properties, \niA\ exceeds tub if i? > is satisfied. Then, when the value of biaxial 
anisotropy is from the vicinity E^, a more rapid thermal variation of rriA results in the N-type dependence 
with one compensation point. As a matter of fact, the total magnetization shows one compensation point 
in which spontaneous magnetization reverses its sign, because at lower temperatures |m^| > ms, while 
at higher temperatures |rn^| < tob (see the curve for E/J = 2.61). In addition, even for more stronger 
biaxial anisotropies the R-type dependences of total magnetization appear {E/J = 2.63). In such a 
case, the total magnetization retrieves its substandard slope from the faster thermal variation of always 
dominating magnetization m^- 

Finally, let us proceed to the discussion of the spin dynamics at non-zero temperatures. The time 
variations of the autocorrelation function C^^j^ are plotted in Fig. 6 for three selected values of biaxial 
anisotropies E/J = 0.1, 0.5 and 2.0. To enable a comparison between the autocorrelation functions at 
various E/J, the relevant temperatures are normalized with respect to their critical temperatures. It 
can be easily realized that the autocorrelation function is not in general a periodic function of time at 
non-zero temperatures. Indeed, C^^j^ arises according to Eq. (28) as a superposition of two harmonic 
oscillations - oscillation with higher angular frequency a;„ = ^ ^/(|)^"+T^)^ another one with lower 
angular frequency uJi = (s)^ (t)^- '^^^ interference between these harmonic oscillations with 
different frequencies and also various amplitudes gives rise to a rather complex time variation of C^^j^, 
which is in general aperiodic, displaying nodes and other typical interference effects (see Fig. 6). The 
dependences drawn in Fig. 6 nicely illustrate also the temperature effect on the spin dynamics. Namely, it 
follows from these dependences, that as the temperature increases, some amplitudes are suppressed, while 
another ones become more robust. Obviously, in the high-temperature region that amplitudes become 
dominant, which coincide to the oscillation with lower angular frequency u)i (see lower panels in Fig. 6). 
In contrast, the amplitudes arising from higher frequency oscillation Wy, dominate at lower temperatures 
(see upper panels in Fig. 6). As far as the influence of biaxial anisotropy is concerned, the stronger the 
ratio E / J , the smaller the difference between both angular frequencies and hence, the more expressive an 
interference effect between them. It is worth mentioning that some particular biaxial anisotropies keeping 
the ratio uju/'^i to be rational, what ensures that the autocorrelation hmction C^^jo is periodic in time 
even at T ^ 0. In any other case, the Cl^^g behaves aperiodically. As a result, in the latter case one 
can impose at best some characteristic time during that the most of spins engaged in the spin dynamics 
change their states, as it apparent from Fig. 6 (al, a2, a3). Although this behaviour is quasi-periodic, 
its characteristic time cannot be confused with the recurrence time r of the former ones, in fact, the spin 
system at biaxial anisotropies giving the irrational ratio oJu/^^l never approaches its initial state again. 

It should be also stressed, that the uniaxial anisotropy D affects the spin dynamics at T ^ 0, as 
well. To illustrate the case, we have depicted in Fig. 7 the time variations of autocorrelation function 
at T/Tc = 1.0 for various uniaxial anisotropies D/J and the ratio E/J = 0.5. Referring to this plot, 
the influence of |0) states on the spin dynamics can be understood more deeply. It is quite evident, 
that a number of the |0) states becomes negligible by taking into account the easy-axis anisotropy (e. 
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g. D/ J = —2.0, Fig. 7a). Really, the quadrupolar moment becomes in this case almost saturated, 
as it can be seen from the time-dependence of C^^^^lt), since C^^jo(O) = q^. Contrary to this, the 
occupation of |0) states becomes crucial when accounting the easy-plane anisotropy (see for instance Fig. 
7d displaying the D/ J = 1.5 case). From the comparison of Fig. 7a-d one can conclude, that positive 
(negative) anisotropy term D reinforces the higher (lower) frequency oscillation Uu {coi) and hence, the 
characteristic time becomes considerably shorter (longer) . Naturally, the observed behaviour results from 
the fact, that the critical temperature gradually falls down as the anisotropy term D/J increases. When 
increasing the ratio D/J, moreover, the oscillation amplitudes are also suppressed, what means, that a 
smaller number of spins changes during the characteristic time their spin states. This striking feature 
clarifies, that the |0) states are not engaged in the spin dynamics so greatly as the | ± 1) states, or even, 
they do not contribute to the spin dynamics at all. 



IV. CONCLUDING REMARKS 



In this article, the exact solution of the mixed spin- 1/2 and spin-1 Ising model on honeycomb lattice 
is presented and discussed in detail. The particular attention has been focused on the effect of uniaxial 
and biaxial crystal- field potentials acting on the spin-1 atoms. As it has been shown, the presence of the 
biaxial anisotropy modifies the magnetic behaviour of studied system in a crucial manner. It turns out 
that already a small amount of the biaxial anisotropy raises a non-trivial spin dynamics and basically 
influences the thermodynamic properties, as well. 

The most interesting finding to emerge here constitutes an exact evidence of the spin tunneling between 
the I ± 1) states in the magnetically ordered phase (OP). Macroscopically, the timneling effect decreases 
the critical temperature for the easy-axis uniaxial anisotropy {D < 0) and also, appreciably depresses 
the magnetization of spin-1 atoms from its saturation value even in the groimd state (see Eq. (26) and 
Fig. 2b). The reduction of critical temperature, as well as magnetization appears apparently due to the 
local quantum fluctuations arising from the biaxial anisotropy. On the other hand, the same quantum 
fluctuations can surprisingly cause an onset of the magnetic long-range order for the extraordinary strong 
easy-plane anisotropies D/J > 1.5. To the best of our knowledge, such a result has not been published 
in the literature before. 

It should be also stressed that there is an interesting correspondence between the model described by 
the Hamiltonian (1) and a similar Ising model with a local transverse magnetic field SI acting on the 
spin-1 atoms only (for a comparison, see Ref.^^): 

3N N 

n = jJ2 ^kf^j + " E ^k- (29) 

However, a similarity between the both Hamiltonians (1) and (29) is not accidental, in fact, when neglect- 
ing the uniaxial crystal-field potential D in the Hamiltonian (1), an effective mapping £^ Q ensures 
the equivalence between (1) and (29). Since this mapping is not related to the magnetic structure in 
any fashion, the appropriate correspondence can be apparently extended to several lattice models. It is 
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therefore valuable to mention, that the magnetic properties of lattice models with the local transverse 
field become a subject matter of many other theoretical works during the last few years^^. Whence, the 
magnetic behaviour of these systems completely resemble that one of their counterparts with the biaxial 
crystal-field potential E only. 

Finally, let us turn back to the origin of biaxial anisotropy. Uprise of this anisotropy term in the 
mixed-spin honeycomb lattice is closely associated with at least a small lattice distortion. To simplify the 
situation, the proposed Hamiltonian (1) accounts the biaxial crystal-field anisotropy, while a diflFerence 
between exchange interactions within various spatial directions of the honeycomb lattice has been, for 
simplicity, omitted. Nevertheless, the developed procedure can be generalized in a rather straightforward 
way also to a model accounting the anisotropic interactions (Ji, J2, J3) within the three non-equivalent 
directions of honeycomb lattice described by the Hamiltonian: 

3N N N 

^ = E '^''^^^(^j + ^ E(^fe)' + ^ E[(4")' - i^in (30) 

(fe,i) fcGB fcGB 

where the nearest-neighbour exchange constant Jkj — Ji, J2 or J3 in dependence on which of the 
three non-equivalent spatial direction it deals. In addition, the biaxial anisotropy strength can be even 
considered as an arbitrary function (linear, quadratic, exponential, logarithmic, ...) of the ratio be- 
tween appropriate interaction parameters: E = /(J2/^i, -/a/Ji)- Therefore, another interesting ques- 
tion arises, namely, whether the Ising model with biaxial crystal-field anisotropy taken as a function 
E = /(J2/J1, J3/J1) can be instable with respect to the spin-Peierls phenomenon. It is quite reasonable 
to assume, however, that under certain conditions the energy gain from the biaxial crystal-field anisotropy 
exceeds the elastic energy related to the lattice deformation and hence, the biaxial anisotropy can lead 
to a spontaneous lattice distortion. To confirm this suggestion, our future work will be directed in this 
way. 
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Figure captions 

Fig. 1 The segment of a mixed-spin honeycomb lattice. The lattice positions of the spin-1/2 (spin-1) 
atoms are schematically designated by full (open) circles, the solid lines label the interactions 
between nearest-neighbouring atoms. The dashed lines represent the effective interaction between 
three outer spin-1/2 atoms, which arise after performing the mapping (9) at fcth site. 

Fig. 2 a) Ground-state phase diagram in the E/ J — D/ J plane; b) Single-site magnetizations (full lines) 
and quadrupolar moment (broken ones) versus the biaxial anisotropy E/J at T = and D/J = 0.0. 

Fig. 3 Time variation of the dynamical autocorrelation function Cf^j^ for various values of biaxial 
anisotropics E/J = 0.1, 0.5, 1.0 and 2.0. Time axis is scaled in h/J units. 

Fig. 4 a) Critical temperature dependence on the uniaxial anisotropy D/J for several values of biaxial 
anisotropics E/J; b) Critical temperature dependence on the biaxial anisotropy E/J for several 
values of uniaxial anisotropics D/J. 

Fig. 5 a) Thermal dependences of the total and sublattice single-site magnetization for D/J = —2.0 and 
E/J = 0.0, 1.0 and 2.0; b) Various temperature dependences of the total magnetization normal- 
ized per one site when the strength of uniaxial anisotropy is fixed [D/J = —2.0) and the biaxial 
anisotropy varies in the vicinity of E'^. 

Fig. 6 Time variations of the dynamical autocorrelation function Cf^j^ when D/J = 0.0 is fixed and 
E/J = 0.1, 0.5 or 2.0. Upper, central and lower panels show the time variation of Cf^j^ at three 
various temperatures, which are normalized with respect to their critical temperatures to ensure 
the same ratio T/Tc = 0.8, 1.0 and 1.1, respectively. 

Fig. 7 The time variations of dynamical autocorrelation function Cf^j^ at critical temperature {T/Tc = 
1.0), E/J = 0.5 and several values D/J = -2.0, 1.0, 1.3 and 1.5. 
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